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A  scheae  for  correcting  for  the  effect  of  ataospheric  refraction  is 
described.  The  scheme  minimizes  the  aaount  of  computation  required  in 
real -tine  applications*  Values  derived  reproduce  observed  data  very  well. 
There  is  a  residual  uncertainty  in  typical  airborne  applications;  using  the 
simplest  scheae  this  amounts  to  about  33/Ur.  A  few  suggested  improvements  to 
the  scheae  reduce  the  uncertainty  below  10  Mr. 
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Recent  years  have  seen  increasing  deployment  of  complex,  airborne  systems 
for  which  highly  accurate  pointing  calculations  are  needed.  These 
calculations  are  used  for  real-time  pointing  of  various  sensors  at  selected 
targets  and  for  reduction  of  observed  target  locations  to  a  common  reference 
frame.  One  of  the  limitations  on  the  accuracy  which  can  be  achieved  in  these 
calculations  is  the  displacement  of  the  apparent  target  location  by 
atmospheric  refraction.  Ground-based  observers  can  use  empirical  data  to 
correct  for  this  displacement.  Airborne  observers  must  predict  the  magnitude 
of  the  effect  for  their  individual  locations  using  some  sort  of  refraction 
model. 

In  a  recent  report  (reference  4)  K-P.  Dunn  describes  one  such  model. 
Dunn's  algorithm  corrects  for  refraction  for  arbitrary  observer  and  target 
locations.  The  calculation  is  geometrically  exact  and  requires  a  rather  great 
amount  of  computation  for  each  correction.  Such  a  large  computational  burden 
may  be  difficult  to  support  in  real-time  applications. 

This  report  discusses  refraction  from  a  somewhat  different  viewpoint  and 
lays  out  a  correction  scheme  specifically  designed  to  minimize  zeal-time 
computation.  In  addition,  the  sources  of  uncertainty  are  evaluated  for 
typical  operating  conditions.  Methods  for  reducing  several  of  the 
contributions  to  the  total  uncertainty  are  proposed. 
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II  ASTRONOMICAL  REFRACTION 


The  precise,  instantaneous  image  displacement  due  to  atmospheric 
refraction  is  very  complicated  and  depends  on  detailed  atmospheric  data  which 
are  not  usually  available,  especially  in  real  time.  As  a  result,  an  exact 
calculation  is  not  possible.  Fortunately,  the  actual  behavior  usually  differs 
only  slightly  from  a  mean  displacement  which  is  simply  an  increase  in  the 
apparent  elevation  of  an  object  with  no  change  in  the  apparent  azimuth.  The 
magnitude  of  this  change  in  elevation  for  an  astronomical  (i.e.,  "infinitely", 
distant)  object  is  called  the  astronomical  refraction.  It  may  be  calculated 
using  various  models  and  may  be  accurately  measured  at  appropriate 
ground-based  observatories.  There  exist  standard  tables  of  astronomical 
refraction;  reference  (1)  is  a  frequently  used  example. 

Astronomical  observations  are  typically  made  at  elevations  greater  than 
about  15°.  In  this  range  the  calculated  refraction  is  insensitive  to  details 
of  the  model  used  and  the  results  are  well  represented  by 

REF  =  R0  cotangent  (elevation)  (1) 

where  RQ  is  a  refraction  parameter  which  depends  upon  atmospheric  conditions 
at  the  observer's  (surface)  location.  At  lower  elevations  the  model  details 
become  important,  and  the  relation  in  EQN  (1)  fails.  In  this  case  recourse  to 
empirical  tables  is  ordinarily  made.  There  are  no  empirical  tables  appropri¬ 
ate  for  airborne  environment  so  a  physical  model  will  be  needed.  The  physics 
needed  for  this  is  straightforward  so  this  approach  is  likely  to  be  reliable, 
but  there  is  a  caveat.  The  empirical  data,  against  which  the  technique  of 
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refraction  modeling  ia  proved ,  consists  of  measurements  made  at  locations 
specifically  chosen  for  their  benign  atmospheric  characteristics.  There 
reeains  soee  chance  that  anoealous  vertical  refraction  (as  in  "mirage” 
phenomena)  and  horizontal  refraction  (azinuth  changes)  nay  be  present. 

Figure  1  shows  the  geometry  of  the  refraction  calculation.  Following 
usual  convention,  a  local  zenith  distance  is  used  instead  of  its  complement, 
the  local  elevation.  Numerical  tables  developed  in  this  report  will,  however, 
be  listed  by  elevation.  An  outward  going  ray*  at  geocentric  distances  r  is 
incident  at  angle  $  on  the  interface  between  media  with  refractive  indices  n 
and  n  +  An.  The  ray  is  refracted  by  Ar  and  proceeds  outwards  to  the  next 
interface  at  r  ♦  Ar  where  it  is  incident  at  angle  ♦  +  A$.  The  flight  over  As 
from  r  to  r  ♦  Ar  subtends  angle  AS  at  the  Barth's  center.  The  physics 
involved  is  Snell's  Law: 

n  sin  ♦  ■  (n  ♦  An)  sin  (♦  +  Ar), 

This  can  be  solved  for  Ar : 

sin  Ar  -  n  sin  ♦  {cos  $  -  [cos2$  *  2  An/n  ♦  (An/n)2]1^2} 

”  n  +  An  (2) 

or  in  the  infinitesimal  limit i 


dr 


dn 

n 


tan$. 


(2a) 


"It  is  somewhat  easier  to  model  the  behavior  of  an  outgoing  ray  instead  of  the 
relevant  Incoming  ray.  Fortunately  Snell's  Law  contains  nothing  about  the 
sense  of  the  propagation  direction  so  the  results  obtained  for  the  outgoing 
ray  will  apply  as  well  to  an  incoming  ray. 
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Addin?  the  interior  angles  of  the  triangle  gives t 


A# 

+  A0  -  Ar  -  0, 

(3) 

d# 

+  d9  -  dT  -  0 

(3a) 

Pinally,  applying  the  cosine  law  to  the  triangle,  leads,  after  considerable 
■anipalation,  to: 

sin  A9  -  {r*  cos**  -  [r'2  cos2*'  -  (r'2  -  r2)]1/2}  (4) 

where 

♦  '  -  ♦  +  A#,  and 
r'  -  r  +  Ar 

This  reduces  to: 

d9  .  tan#  (4a) 

Inserting  (4a)  and  (2a)  into  (3a)  leads  to 

d(nr  sin  #)  ■  0. 

Evaluating  the  product  at  the  observer's  location  gives  the  useful 
relationship: 

nr  sin  #  »  n0ro  sin  #Q  (5) 

Note  that  #Q  is  the  observed  (apparent)  zenith  distance  of  the  ray  arriving 
at  the  observer's  location. 

Although  equations  (2),  (3),  and  (4)  are  the  exact  equations  for  finite 
steps,  and  equations  (2a),  (3a),  and  (4a)  give  only  approximate  values,  it 
should  be  noted  that  the  finite  steps  are  already  an  approximation  of  the 
continuous  atmosphere.  In  addition,  equations  (2),  (3),  and  (4)  comprise  a 
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set  of  coupled  transcendental  equations  —  a  nasty  business.  Integration  was 

therefore  carried  out  using  the  infinitesinal  forms.  These  were  recast  to 

avoid  problems  with  tan  $  at  small  elevation  angles  ($  ■  */2) . 

Defining  a  refraction  parameter 

„  _  1  dn 

”  T?  dr 

The  equations  for  the  model  become: 

Ar  m  y  sin$  As, 

A$  -  (Y  -  l/r)  sin$  As,  and 
Ar  *  cos$  As. 

These  equations  were  then  used  to  integrate  outward  from  the  observer  (at 
rQ,  <PQ)  to  a  point  at  which  negligible  additional  refraction  is  produced. 
Initially  the  step  3ize  As  was  set  very  small  and  the  upper  limit  on  r  very 
large.  The  step  was  then  increased  and  the  upper  limit  decreased  until  an 
error  of  about  1  Ur  was  produced.  This  point  was  arbitrarily  chosen  as  the 
optimum  point  in  the  trade-off  of  accuracy  vs.  speed. 

The  "natural"  input  parameter  for  calculation  of  T  is  the  apparent  zenith 
distance  $Q.  When  treated  as  a  function  of  the  apparent  zenith  distance, 
the  refraction  will  be  denoted  rj.  By  iterating  on  the  input  angle  it  is  also 
possible  to  produce  values  of  T  treated  as  a  function  of  ($0  +  T)  —  i.e., 
the  true  zenith  distance  of  an  astronomical  object.  When  this  is  done,  the 
refraction  will  be  denoted  Each  version  has  its  place.  Throughout  most 

of  this  report  T  will  be  written  without  subscript  because  the  choice  is 
irrelevant  to  the  point  under  discussion.  In  fact,  the  numerical  examples 
were  obtained  with  version  I*i,  except  for  a  single  ?2  value  on  page  24  . 
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In  order  to  calculate  the  refraction  it  is  necessary  to  derive  values  of 
the  refraction  parameter  Y.  Several  possibilities  exist.  The  simplest  would 
be  to  invent  a  model  of  n  as  a  function  of  r.  A  better  approach  would  be  to 
model  the  atmospheric  pressure  and  temperature  profiles  and  calculate  n  from 
these.  The  approach  taken  here  was  to  use  the  temperature  and  pressure 
profiles  given  in  the  set  of  standard  supplemental  atmospheres  (reference  2) 
and  calculate  n  from  these.  The  supplemental  atmospheres  represent  a 
combination  of  physical  theory  and  empirical  data;  they  are  given  for  several 
latitudes  and  for  summer  and  winter  conditions.  The  refractive  index  was  then 
calculated  using  currently  adopted  relationships  (reference  3): 

(n-1  )  x  106  -  272.6  q(A)  f(p,t) 

where 

q(A)  -  0.236  ♦  108.2/(146  -  I/*2)  ♦  0.937/(41  -  1/A2) 
with  A  in  Urn,  and 

f(  m  p  [l  +  p( 1.049  -  .01 57t)  x  10“6] 

P'  720.893  (1 .  +  ,00366t) 

with  p  in  nun  Hg  and  t  in  “C.  The  wavelength  factor  goes  to  unity  for 
increasing  wavelength;  q(.55)  »  1.019;  q(10)  -1  <<  .001.  The  thermodynamic 
state  factor  is  unity  at  760  mm  Hg  and  15°C.  Seven  sets  of  n  have  been 
compiled:  one  for  the  year-round  tropical  atmosphere  (TRSW);  two  for 
midlatitudes,  a  summer  (MLSU)  and  a  winter  (MLWI);  and  four  for  subarctica,  a 
summer  (SASU)  a  mean  winter  (SAWI),  a  cold  winter  (SAWC),  i.-.d  a  warm  winter 
( SAWW ) .  The  index  of  refraction  was  calculated  for  1  km  steps  up  to  a  height 
of  50  km.  In  the  ray  tracing  calculation  a  four-point  interpolation  was  used 
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to  give  n  and  its  derivative  as  functions  of  h.  Beyond  50  km  a  constant 
exponential  decrease,  which  was  derived  from  the  40-50  km  points,  was  assumed. 

Zt  would  be  highly  desirable  to  test  the  refraction  model  described  above 
against  empirical  data  obtained  under  a  variety  of  conditions.  Onfortunately 
the  only  such  data  available  is  the  accumulated  ground-based  experience, 
obtained  under  a  very  narrow  range  of  conditions.  Table  Z  shows  the  results 
of  refraction  calculated  using  the  MLSD  and  MLWZ  atmospheres  with  the  observer 
at  sea- level.  Zn  order  to  facilitate  comparison  with  published  discussion  of 
refraction  and  measured  values  of  refraction,  the  astronomers'  angle  unit,  the 
second  of  arc,  is  used;  1*  ■  4.85  Ur.  Also  listed  is  the  empirical  data  given 
in  reference  (1).  The  agreement  is  satisfactory.  The  last  column  shows  the 
refraction  calculated  using  Bqn  ( 1 ) .  Good  agreement  is  obtained  above  40°  but 
serious  discrepancies  have  developed  by  10°  elevation. 
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TABLE  Z 


ASTRONOMICAL  REFRACTION  AT  SEA-LEVEL 


A Elevation  (arc- sec) 


Apparent 

Elevation 

Calculated« 

Winter 

Mid-Latitude 

Suaaer 

Allen 

AQ 

Ro* 

Cot(e) 

45* 

61 

56 

59 

59 

40 

72 

67 

70 

71 

35 

87 

80 

84 

85 

30 

105 

97 

101 

103 

25 

130 

120 

125 

127 

20 

166 

153 

159 

163 

15 

224 

206 

215 

222 

10 

334 

307 

319 

337 

8 

412 

378 

394 

423 

6 

533 

488 

509 

565 

4 

741 

674 

706 

849 

3 

911 

824 

865 

1133 

2 

1163 

1043 

1103 

1701 

1 

1566 

1384 

1481 

3403 

0.5 

1874 

1638 

1760 

6945 

0.0 

2279 

1963 

2123 
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Linearity  and  Separability 

Apart  from  the  simple  geometric  functionality,  EQN  (1)  carries  two 
implications  that  are  ordinarily  assumed  true  even  if  a  more  complicated 
elevation  dependence  is  used.  These  are  that  the  refraction  depends  linearly 
on  some  refraction  parameter  (RQ,  which  may  be  an  arbitrary  function  of 
temperature,  wavelength,  etc.)  and  that  the  geometrical  and  refractive  effects 
are  separable.  That  this  cannot  be  strictly  true  is  easily  seen  by  examining 
the  details  of  refraction  in  the  atmosphere.  Consider  two  rays  outward  bound 
at  0,  with  the  same  zenith  distance  but  having  different  wavelengths  and 
consequently  seeing  different  refractive  index.  At  the  first  interface  the 
ray  seeing  the  higher  refractive  index  will  experience  greater  refraction  in 
proportion  to  its  greater  refractive  index.  At  the  second  interface  this  same 
ray  will  experience  greater  refraction  both  because  of  its  higher  index  and 
also  because  of  its  greater  zenith  distance.  The  total  refraction  must  thus 
increase  faster  than  linearly  with  Rq.  Because  the  geometrical  dependence 
is  not  linear  with  elevation,  the  magnitude  of  the  supralinearity  must  depend 
on  the  elevation,  and  so  the  effects  couple. 

This  behavior  can  be  seen  by  exercising  the  model.  Using  MLWI,  a 
ground-based  observer,  elevations  of  5°  and  10°,  and  values  of  q(X)  of  1.0  and 
1.2  (physically  implausible)  both  effects  are  produced.  At  5°  the  total 
refraction  for  q  >  1.2  is  20.3%  higher  than  for  q  *  1 .0  instead  of  the  linear 
20%.  This  implies  generation  of  an  error  of  9.4  Vr  if  the  linear  assumption 
is  used.  At  10°  the  increase  has  dropped  to  20.1%.  In  addition  to  verifying 
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the  coupled,  nonlinear  behavior,  these  numbers  show  that  the  effects  are 
small.  Ne  shall  therefore  adopt  the  fora  of  EQN  (1),  using  the  model  to  get 
the  geometrical  relationship  at  a  nominal  value  of  Rq  and  then  making  a 
linear  adjustment  for  actual  wavelength  and  environmental  conditions. 

Use  of  Tables  and  Interpolation 

Although  the  refraction  calculation  was  adjusted  to  run  as  quickly  as 
possible  without  producing  significant  error,  it  still  takes  about  10 
milliseconds  to  run  on  the  VAX  11/750.  This  is  a  burden  in  real-time 
operation  and  should  be  avoided  if  possible.  An  alternative  is  use  of  a  table 
and  some  more  or  less  complicated  interpolation  scheme.  The  simplest  case 
would  be  regularly  spaced  sample  points  with  a  linear  interpolation.  Table  II 
shows  the  results  of  such  a  scheme.  First  the  model  calculation  was  used  to 
give  refraction  values  at  half -degree  intervals.  The  SASU  atmosphere  at  11.3 
km  was  used  in  the  calculation.  Next  values  for  half-odd-integer  values  of 
elevation  were  obtained  by  averaging  the  surrounding  integer  elevation 
values.  The  resulting  error  is  given  in  the  last  column. 

Because  of  the  nature  of  the  dependence  of  refraction  on  elevation, 
linear  interpolation  always  overestimates  refraction.  The  errors  at  the  mid¬ 
points  can  thus  be  reduced  by  subtracting  half  the  error  from  the  values  for 
whole  degree  elevations.  Table  III  shows  such  a  revised  table,  interpolated 
values,  model  calculations,  and  interpolation  errors.  As  expected,  the  error 
at  the  half  integer  elevation  is  reduced  at  the  expense  of  adding  an  error  at 
the  integer  elevations. 
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TABLE  II 

LINEAR  INTERPOLATION  ERROR 


Apparent 

Calculated 

Interpolated 

Interpolation 

Elevation 

Refraction 

Refraction 

Error 

(deg) 

( . 

- Wr - - 

- - ) 

9.0 

461 .7 

8.5 

487.2 

488.6 

+1 .4 

8.0 

515.5 

7.5 

547.1 

549.1 

+2.0 

7.0 

582.6 

6.5 

622.7 

625.5 

+2.8 

6.0 

668.4 

5.5 

720.7 

724.9 

+4.2 

5.0 

781 .3 

4.5 

852.0 

858.4 

+6.4 

4.0 

935.5 

3.5 

1035.5 

1045.9 

+10.6 

3.0 

1156.2 
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TABLE  XZZ 


INTERPOLATION  ERROR  WITH  ADJUSTED  TABLE 


Apparent 

Elevation 

(deg) 

Adjusted 
Refraction 
( . 

Interpolated 

Refraction 

Calculated 

Refraction 

Interpolation 

Error 

. ) 

4.2 

900.9 

900.3 

+0.6 

4.1 

916.1 

917.6 

-1 .5 

4.0 

931.4 

931.4 

935.5 

-4.1 

3.9 

953.2 

954.0 

-0.8 

3.8 

975.0 

973.2 

+1 .8 

3.7 

996.8 

993.1 

+3.7 

3.6 

1018.6 

1013.8 

+4.8 

3.5 

1040.4 

1035.3 

+5.1 

3.4 

1062.1 

1057.6 

+4.5 

3.3 

1083.9 

1080.8 

+3.1 

3.2 

1105.7 

1104.9 

+0.8 

3.1 

1127.5 

1130.0 

-2.5 

3.0 

1149.3 

1149.3 

1156.2 

-6.9 

2.9 

1182.3 

1183.5 

-1  .2 

2.8 

1215.3 

1211  .9 

-3.4 

14 


Although  the  Interpolation  error  la  strictly  deterministic  in  nature  we 
may  get  some  feeling  for  its  importance  by  analysing  it  as  a  random  error. 
Following  the  procedure  above,  the  RMS  interpolation  error  in  the  range  0°  to 
1 •  elevation  is  30  Ur  or  about  1.2%  of  the  refraction  at  0.5°  elevation, 
this  is,  at  worst,  not  the  dominant  source  of  error  and  is,,  in  any  event, 
negligible  above  a  few  degrees  of  elevation.  Thus  use  of  a  linear 
interpolation  on  a  table  of  values  at  1*  intervals  provides  satisfactory 
accuracy  and  vastly  improved  computational  speed.  Construction  of  appropriate 
tables  should  be  carried  out  in  the  mission  planning  phase. 

Sensitivity  to  observer  Height 

The  height  h  of  the  observer  is  an  input  parameter  for  the  Refraction 
model.  As  in  the  case  of  the  refraction  parameter,  there  is  no  reason  to 
expect  the  exact  dependency  of  total  refraction  upon  h  to  be  simple.  For 
small  enough  departures  from  nominal  height  hg  we  can,  however,  expect  to  be 
able  to  use  an  equation  of  the  form 

T (h,e)  -  f (h,e)  T  (hQ,e) 

in  order  to  simplify  calculation  of  T  when  a  range  of  h  is  anticipated.  In 
fact,  for  typical  airborne  operational  altitudes  good  results  can  be  obtained 
with  a  very  simple  expression  for  f: 

f(h,e)  -  fQ  +  f ' (h  -hQ) . 


L 
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Using  SASCJ  and  Kq  -  11.3  km  the  beet  fit  in  the  range  10.3  <  h  <  12.3  km  is 
given  by 

f  -  1.0035455 
o 

f  -  -0.1494091 

For  e  *  10°  the  "RMS  error”  introduced  by  this  aethod  ia  about  1.3  tir.  The 
typical -uncertainty  in  the  observer's  height,  about  30  a,  propagates  through 
either  the  full  calculation  or  the  linearized  correction  to  give  an 
uncertainty  of  1.9  Ur  at  the  saae  nominal  point. 

Variation  with  Latitude 

Reference  (2)  gives  model  atmospheres  for  latitude  15”,  30”,  45”,  60”, 
and  75°.  These  differ  not  only  in  the  mean  temperature,  but  also  in  subtler 
waye  such  as  scale  height.  In  order  to  determine  the  effect  of  latitude  on 
the  refraction,  calculations  were  carried  out  for  tropic,  mid-latitude,  and 
sub-arctic  atmospheres.  Table  IV  shows  the  results  for  the  tropic  atmosphere 
and  for  the  middle  of  the  ranges  spanned  by  the  several  mid-latitude  and 
sub-arctic  atmospheres.  The  calculations  were  carried  out  for  q  -  1 .0  and 
h  *  11.3  km.  The  differences,  while  not  large,  are  significant  at  low 
elevations.  Clearly  it  is  worth  the  trouble  to  chose  the  appropriate 
atmosphere.  For  intermediate  latitudes  there  is  no  choice  but  to 
interpolate.  This  is  clearly  not  physically  correct,  but  seems  unlikely  to 
introduce  more  than  a  few  Ur  error  at  e  »  10°. 
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TABLE  XV 


LATITUDE  EFFECT 


Apparent 

Elevation 

TR 

Refraction  (p r) 

ML 

SA 

20° 

220 

211 

197 

18 

246 

235 

220 

16 

278 

266 

248 

14 

318 

305 

285 

12 

371 

355 

333 

10 

443 

424 

398 

8 

547 

523 

491 

6 

708 

679 

638 

4 

988 

950 

893 

2 

1556 

1510 

1426 

0 

2966 

2997 

2924 
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Variation  with  Saaaon 

Reference  (2)  gives  four  sub-arctic  atmospheres,  one  for  suaaer,  and 
three  for  winter.  The  accoapanying  text  points  out  that  the  winter-warm  and 
winter-cold  ataosphares  do  not  represent  extremes  of  a  continuous  population 
but  rather  are  typical  naabers  of  an  essentially  bimodal  distribution.  Table 
V  shows  the  refraction  for  the  four  atmospheres;  again  q  *  1 .0  and  h  ■  11.3ka. 
The  seasonal  effect  can  be  seen  to  be  comparable  to  the  latitude  effect,  so 
again  it  is  worthwhile  to  use  the  appropriate  ataosphere  in  the  calculations. 
In  the  case  of  intermediate  seasons,  the  biaodal  distribution  in  winter  makes 
hash  of  the  idea  of  interpolation.  A  plausible  course  might  be  to  use  SASU  in 
suaaer,  choose  somehow  (weather  data,  perhaps)  between  SAWW  and  SAWC  in 
winter,  and  use  a  table  aidway  between  SASU  and  SAWI  the  rest  of  the  year. 

This  will  undoubtedly  introduce  soae  error.  Judging  from  the  table,  this 
error  should  be  negligible  above  20*,  be  about  4-5  ur  at  10*,  and  grow  to 
about  30  Ur  by  0*. 

Effect  of  Weather 

An  airborne  measurement  of  an  atmospheric  parameter,  say  temperature, 
will,  in  general,  differ  from  the  corresponding  value  in  the  atmospheric 
profile  used  in  the  refraction  calculation.  For  example,  consider  a  winter 
flight  at  11.3  km  with  the  SAWC  atmosphere  chosen  as  appropriate.  At  that 
height  the  model  profiles  predict  about  225K,  222K,  21 7K,  and  21 7K  for  SU,  WW, 
Wl,  and  WC  respectively.  A  measurement  of  219K  would  have  contributions  from 
measurement  error,  a  local  perturbation  from  the  nominal  profile,  and  a  global 
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TABLE  V 


SEASONAL  VARIATION  OF  REFRACTION 


Refraction  (Ur) 


Apparent 


Elevation 

SASU 

SAWI 

SAWN 

SAWC 

20° 

206 

192 

188 

192 

18° 

230 

214 

210 

214 

16* 

260 

242 

237 

242 

14° 

298 

278 

272 

278 

12* 

348 

324 

318 

324 

10° 

41  o 

387 

380 

387 

8° 

514 

479 

469 

479 

6" 

666 

622 

609 

622 

4» 

933 

872 

853 

872 

2° 

1489 

1396 

1363 

1394 

0 • 

3034 

2876 

2813 

2860 
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shift  of  the  entire  temperature  profile.  These  are  completely  different 
physical  situations;  how  should  the  +2K  temperature  shift  be  apportioned? 

Even  if  the  first  two  contributions  could  be  ruled  out,  the  situation  is  no 
better.  Should  the  refraction  be  increased  by  interpolation  between  WC  and  SO 
(217  and  225)  or  decreased  by  interpolation  between  WC  and  WW  (217  and  222)? 

Zt  appears  that  it  will  not  be  useful  to  attempt  to  correct  for  weather.  This 
night  introduce  an  uncertainty  of  about  the  same  magnitude  as  the  uncorrected 
seasonal  variation.  More  likely  the  effect  will  be  considerably  smaller 
because  airborne  operations  take  place  near  or  above  the  tropopause. 


IV  THE  POSITION  CORRECTION 


The  astronomical  refraction  T  and  the  apparent  change  in  elevation  e  are 
equal  only  for  stars.  For  nearby  objects  the  ray  which  will  reach  the 
observer  after  refraction  does  not  start  out  parallel  to  the  radius  vector 
from  object  to  observer.  The  geometry  is  shown  in  Figure  2.  A  ray  (dashed 
line)  starts  from  the  target  point  T  towards  the  apparent  observer  O',  with 
local  zenith  distance  $ .  The  ray  is  refracted  through  angle  T  and  arrives  at 
the  observer  O  with  apparent  zenith  distance  The  true  (geometric) 

zenith  distance  of  the  target  is  C.  If  enough  of  the  geometry  is  known,  the 
elevation  correction  e  may  be  calculated.  Two  cases  occur:  (i)  an  object  is 
observed  at  <f>0  and  it  is  desired  to  find  the  true  zenith  distance  C  (or 
elevation,  * /2  -  C)>  in  this  case  the  appropriate  argument  for  T  is  clearly 
the  apparent  zenith  distance,  which  is  known,  (ii)  the  true  zenith  distance  of 
an  object  is  calculated,  and  it  is  necessary  to  determine  the  apparent 
position  at  which  to  point  the  telescope.  In  the  second  case  the  argument  for 
T  is  a  problem.  The  value  of  $Q  is  the  unknown,  and  $0  4  T  -  C  only  for 
very  large  distances.  As  a  result  it  is  necessary  either  to  iterate  or  to 
accept  an  approximation. 

Case  (i)  is  appropriate  for  the  reduction  of  observed  locations  to  a 
standard  frame.  Angle  COT  can  be  seen  to  be  w  -  -  e,  angle  OTC  is  $  -  T 

+  e,  and  consequently  angle  OCT  is  -  $  +  T .  The  observed  zenith  distance 
is  known  and  r  can'  be  found  as  a  function  of  The  trigonometric 
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equations  plus  equation  (5)  fora  a  set  which  can  be  solved  for  £.  After 
significant  manipulation: 

n  cos  r  -  n  cot  +  sin  I*  -  n 
e  m  o  _ o  _ _ 

n  cot  4  -  n  cot  $  cos  T  -  n  sin  t 

o  o  o 


This  can  be  recast  in  a  somewhat  more  illuminating  fora  by  eliminating  $  in 
favor  of  p,  the  geometric  range  to  the  target,  and  solving  for  the  difference 
between  the  refraction  and  the  elevation  correction.  Assuming  the  refractive 
index  at  the  target  is  unity: 

r 

sin(T  -  £)  -  -jp.  [nQ  sin*0  -  sin(+Q  +  I*)] 


Note  that  this  agrees  with  the  astronomical  limit:  for  p  *  e  ♦  I*. 

Case  (ii)  is  appropriate  for  calculation  of  an  apparent  position  from 
physical  data.  In  this  case  the  triangle  OCT  is  completely  known.  Combining 
the  trigonometric  equations  with  equation  (5)  leads  to: 

nr  sin  (X  +  T)  -nr  sin  C 
o  o 

tan  6  m  . ;  \  f . .  _ 

nr  cos  (A  +  r)  -  n  r  cos  C 
o  o 


Again,  a  handier  relationship  may  be  obtained  for  T  -  e.  For  n  ♦  1 , 


tan(r-e) 


n  sin  (C  -  T)  -  sin  C 
_ o _ 

cos  (C  -  1*)  +  cos  C  +  p/r  ' 


or 


(6) 


[  n  sin  (S  -  T )  -  sin  c] 


(f  -  E) 


(6a 


As  was  noted  above,  neither  of  the  possible  arguments  for  T  is  known.  A 
nuaerical  exaaple  shows  how  to  deal  with  this  and  also  puts  the  entire  problem 
in  proper  perspective.  For  input,  use  the  following: 

Observer  Height:  h  -  11.3  ka 
Atmosphere:  SAWI 
True  Zenith  Distance:  C  -  80s 
Range  to  Target:  P  «  400  ka 

Wavelength  Factor:  q  ■  1 .0 

Observer  Geocentric  Distance:  r0  ■  6375.15  ka 

The  first  approximation  for  the  refraction  can  be  obtained  using  C  as 
the  argument  either  in  ri($0)  or  I*2(^Q  +  I*).  These  give: 

r i ( C )  >  389.2  Ur  and 
1*2  <C ) .  -  388.4  Mr 

TJie  remainder  of  the  example  will  continue  with  the  Tj  version.  Next, 

T  -  c  can  be  calculated  either  from  EQN  (6)  or  EQN  (6a): 

r  -Cm  26.8  Ur  (EQN  6) 

T  -  E  -  26.6  Ur  (EQN  6a) 
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T-_v 


From  this  a  new  approximation  to  $0  ■  C  -  e  may  be  obtained: 

♦0  -  79.979236", 

leading  to  a  new  T 

Ti  -  388.4  Ur, 

and  a  new  (T  —  c ) 

T  -  c  -  26.7  Ur, 

a  new  40 

♦0  -  79.979276", 

a  new  T 

T!  -  388.4  Ur, 

and  so  on.  Clearly  this  is  far  enough.  The  following  conclusions  are  now 
apparent: 

-  The  first  approximation  using  V2  or  the  second  approximation  using  Tj 
are  both  adequate. 

-  Equation  (6a)  is  sufficiently  accurate  (unless  p  is  very  small) 

-  The  elevation  correction  e  is  approximately  equal  to  the  refraction  T . 

This  last  conclusion  is  strengthened  when  greater  ranges  are  considered.  For 
P  =t  600  km,  (T  -  e)  >  17.8  Ur,  and  for  P  -  800  km,  (I*  -  e)  «  13.4  Ur. 

Airborne  systems  are  not  ordinarily  controlled  in  azimuth  a  and  zenith 
distance  z,  but  in  an  inertial  coordinate  system.  For  this  memo  the 
coordinate  system  will  be  taken  to  be  the  astronomical  right  ascension  (a), 
declination  (5)  system,  although  there  will  probably  be  a  flight-constant 
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offset  froa  true  <*.  Thus  it  trill  be  necessary  to  calculate  corrections  to  a 
and  5  which  are  the  result  of  the  correction  e  to  the  zenith  distance.  The 
relation  between  the  two  coordinate  systems  involves  both  a  constant  tilt 
which  depends  on  the  astronomical  latitude  B  and  longitude  A,  and  a  slip  which 
involves  the  Greenwich  sidereal  time  9g.  An  intermediate  coordinate  system 
-  hour  angle  h,  declination  6  is  usually  defined.  This  system  has  the  same 
pole  as  the  (a, 6)  -  system  but  uses  the  observer's  meridian  as  the  zero  of  the 
cyclic  coordinate  h,  which  is  a  left-handed  angle.  The  relationships  are: 

h  ■  8g  +  A  -  «,  and 
6  -  « 

Figure  3  shows  the  relationship  between  (h,6),  and  (a,z).  In  this  figure 
Z  is  the  observer's  zenith,  P  is  the  north  celestial  pole,  and  T  is  the  two- 
dimensional  location  of  the  target.  The  systems  are  related  by  a  fixed 
rotation : 

sin  z  cos  a\  /-sin  B  0  cos  b\  /  cos  6  cos  h\ 

sin  z  sin  a  I  =  I  0  -t  0  I  cos  A  sin  h  j 

^  cos  z  j  \  cos  B  0  sin  B  J  1  sin  5  I 

which  can  be  solved  to  give  (a,z)  as  a  function  of  (h,6)  and  B  or  (h,6)  as  a 

function  of  (a,z)  and  B.  Alternatively,  the  spherical  triangle  PZT  may  be 

solved  to  provide  these  equations. 

The  straightforward  way  to  calculate  the  corrections  to  h  and  6  is  to 
solve  for  (a,z)  corresponding  to  (h,6),  correct  for  refraction  by  a'  ■  a  and 
z*  =*  z  -  £ ,  and  finally  solve  for  (h',5')  corresponding  to  (a',z').  This  is 
clearly  an  inappropriately  large  amount  of  computation.  It  is  possible  to 
reduce  the  amount  of  computation  somewhat  by  taking  advantage  of  the  small 
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Figure  3.  Relationship  between  the  celestial  and  horizon  coordinate  systems 


angle  approximations  for  functions  of  G .  This  is  furthsr  facilitated  by  use 
of  an  auxiliary  quantity,  the  parallactic  angle  p  in  Figure  3.  The 
parallactic  angle  nay  be  found  using  the  following  sequence  of  calculations: 

cos  z  -  sin  B  sin  6  +  cos  B  cos  i  cos  h 
sin  x  m  +[i  -  cos2z]V2 
cos  p  -  (sin  B  -  cos  z  sin  6)/sin  z  cos  & 
sin  p  •  cos  B  sin  h/sin  z 

then 

6*  -  6  *  e  cos  p,  and 

o'  -  a  -  h  -  h*  (7) 

“  e  sin  p/cos  5* 

At  z  *  80*,  a  -  300*,  and  G  -  400  Ur,  use  of  these  approximations  introduces 
and  error  of  about  a  hundreth  of  a  microradian.  Further  simplifications  may 
be  appropriate  in  specific  applications.  These  include  use  of  encoder  values 
for  z  Instead  of  their  calculation,  use  of  mean  or  constant  values  for  B  and 
p,  etc.  Evaluation  of  the  error  introduced  by  such  measures  is  beyond  the 
scope  of  this  study. 
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V  SUMMARY 

The  following  scheme  Is  suggested  for  correcting  calculated  pointing 
angles  for  refraction  in  an  operational  environment. 

*  During  mission  planning,  calculate  tables  of  astronomical  refraction  as  a 
function  of  true  elevation  (i.e.,  r^) .  The  input  parameters  should  span 

the  anticipated  ranges  of  observer  height,  latitude,  and  season.  The  tables 

should  be  adjusted  to  correct  for  the  effect  of  linear  interpolation  (see  p. 

12).  Values  of  n  should  also  be  calculated, 
o 

*  During  mission  planning,  produce  a  single  refraction  table  and  a  single 
value  of  nQ  from  the  results  of  the  previous  step.  These  should  be 
evaluated  for  the  nominal  observer  height,  and  "interpolated"  for  specific 
latitude  and  time  of  year.  Values  of  dnQ/dh,  fQ,  and  f'  (see  p.  16) 
should  also  be  obtained  from  the  variations  over  the  range  of  observer 
heights . 

•  In  near-real  time,  obtain  the  appropriate  values  of  for  the 
relevant  elevation  by  linear  interpolation  from  the  table.  Correct  r and 
nQ  for  the  actual  observer  height. 

•  In  near-real  time,  calculate  the  elevation  correction  e  using  EQN  (6a) 
and  the  pointing  corrections  Act  s  A5  using  EQN  (7). 
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Correction  for  atmospheric  refraction  on  the  basis  of  this  scheme 
involves  several  assumptions  and  approximations.  Table  VI  summarizes  these; 
magnitudes  of  uncertainties  are  evaluated  for  10s  elevation,  and  observer 
height  11.3  km.  Uncertainties  in  the  refraction  correction  all  have  very  low 
frequency  spectra.  The  first  entry  is  a  reminder  that  there  is  no  evidence 
which  allows  us  to  rule  out  the  possibility  of  anomalous  refraction.  The  most 
problematical  of  the  remaining  uncertainties  is  the  weather  effect.  The  value 
is  in  fact  little  better  than  a  guess  at  the  upper  limit.  If  the 
uncertainties  are  combined  (RSS)  it  is  seen  that  even  a  crude  treatment  ia 
satisfactory  for  many  airborne  system  missions,  giving  a  total  uncertainty  of 
about  33Mr.  In  this  case  the  weather  effect  is  just  one  of  several 
contributors,  so  reduction  of  this  figure  is  of  little  consequence.  It  is 
probably  desirable  to  improve  this,  or  equivalently  to  produce  this  acceptable 
error  at  lower  elevations.  This  ia  particularly  true  in  light  of  the 
relatively  low  burden  imposed  by  the  improved  calculations.  In  that  case,  the 
size  of  the  weather  induced  uncertainty  is  crucial.  If,  as  was  suggested  in 
Section  III,  the  figure  listed  is  a  great  overestimate  of  the  actual  effect  of 
weather,  then  the  total  error  can  be  reduced  to  less  than  10  Mr. 
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TABLE  VI 


SUMMARY  OF  UNCERTAINTIES  AND  ERRORS 
IN  THE  REFRACTION  CORRECTION 


Source 


Assumption  that  refractive  index  is  monotonic 
function  of  height  only 


Assumption  of  linearity  and  separability 


Use  of  table  and  interpolation 


Sensitivity  to  observer  height 


Latitude 

Uncorrected 

Effect 

Interpolated 

Seasonal 

Uncorrected 

Variation 

Interpolated 

Heather 

Calculation 

Set  equal  to  T 

of  e 

Equation  ( 6a ) 

Calculation  of  A a,  A5 


Magnitude 

Unknown 

<  1  Ur 

<  1  Ur 

»  2.5  ur 

“  20  ur 

*  2  ur 

-  15  Ur 

•  5  ur 

<  15  Ur? 

-  15  ur 

<  1  Ur 

«  1  Ur 


2)  A.  E.  Cole,  A.  Court,  and  A.  J.  Kan  tor,  "Model  Atmospheres",  in  Handbook 
96  Geophysics  and  Space  Environments,  (AF  Cambridge  Research  Laboratories, 
1965). 


3)  D.  E.  Gray,  ED.,  American  institute  of 
New  York,  1972). 


lies  Handbook,  (McGraw  Hill, 


4)  K-P.  Dunn,  "Atmospheric  Refraction  Error  and  Its  Compensation  for  Passive 
Optical  Sensors",  Technical  Report  686,  Lincoln  Laboratory,  M.I.T. , 

4  June  1984). 
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simplest  scheme  this  amounts  to  about  33  pr.  A  few  suggested  improvements  to  the  scheme  reduce 
the  uncertainty  below  10  pr. 
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